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Abstract 

We write a space-time Feynman Path Integral representation for scattered elastic wave 
fields from a weakly anisotropic non-homogeneity. 



1 Introduction 



It has been an useful tool to study waves of various physical nature, the standard by now oper- 
atorial method of Feynman path integral ([!]). However the generahzation of the path integral 
method to wave propagation in solids as well as liquids needs proposals for extension of the path 
integral technique to a full elastic-anisotropic medium, a subject of the utmost importance to 
Geophysics (see ref. [2] for insights into the origin of the method of path summation for elastic 
waves). We aim in this note to contribute to this still unsolved problem of anisotropic wave 
propagations by considering a sort of approximate - first order perturbation Feynman path 
integral representation for wave scattering through a weakly anisotropic scattering potential 
(second reference in [1]). 



2 The Anisotropic Path Integral Representation 

Let us start our analysis by considering a wave equation associated to an "acoustic wave" 
scattered by an anisotropic center in with a scalar space-time dependent damping 

-Q^U\x,t) = Ci^k,{x)-^-^^U\x,t) - u{x, t)-U'{x,t) + 5{x - x')5{t - t')ei (1) 

where the time range is the whole real line — oo < i < oo, Cj are the standard canonical 
vectorial basis of R^, W{x,t) is the vector position of the wave field under scattering and 
lijke{x) is a spatially-dependent anisotropic medium tensor carrying the anisotropic scattering 
center properties {'fijkiix) = ^jiki{x) = 'yijik{x) = ^ktij{x)) 

Cijkt{x) = [CQ{x)SikSji + €'jijke{x)]. (2) 

Note that Cq{x) is the acoustic unperturbed wave vecolicity scalar function and ^ijke{x) is 
a fourth-order "elasticity" tensor of the scatering anisotropic center and £ is a small parameter 
associated to the intensity of this wave field weakly perturbation. The positive definite function 
^{x, t) describe our "losses" due to the existence of damping. 

The Green function associated to our coeficients variable anisotropic wave equation Eq.(l) 
is supposed to have the operatorial-distributional rule above written as its definition (second 



reference in [1] and the convention 'jajbev-'v^ = jUv^v^I) in quantum mechanical Dirac bra-ket 
notatioij^ 
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Since e is small and the "scalar acoustic tensor" Cg (^)'^ife^i^ f^Hy diagonal, one can proceed 
as in refs ([1] and [2]) and arrive at the full space-time wave field propagation generated by a 
source at the point x' and at the time t'. (See Appendix B for the calculational details) 
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where we have introduced the Anisotropic Polarizations vector Factor directing the wave prop- 
agation along the anisotropics axis 
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Note the use of an proper-time ordered - matrix indexes {p, q) on the above written anisotropic 
factor. 



^We have used the L. Schwartz distributional integral representation for a > 0: 
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That is our general result. However, altliouglit a somewhat intrincate object to be analiti- 
cally evaluated ([1]), in practical lossless geophysical/signal detection-eigeneering applications 
(with iy{x,t) = 0), one is always interested in the so called "short-time arrival signals". 

At this WKB like approach for our path integralj^, Brownian like path trajectories {r{cr)} 
roughs/ nondifferentiable curves are dominated by geometrical rays paths, i.e., those classical 
different iable paths are dominant paths contributing to our path integral Feynman measures 
(including the "classical" weights!). We have thus the standard results ([l]-first reference, 
chapter 24) 

D^[t(a)]e(i/o(^)^) = ^^;= (6) 
t{o)=t' VAnis 

f(o)=x' (y/Inis)^ 



and 



i^ik - Sik) ~ -is 



where the geometrical ray path satisfies the calssical ODE's rays equation of classical optics 
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Cl{r{a)) 

and the geodesic lenght connecting the points of observation x and x' is given for a small 
desviation condition by the quadratic object 



A2(x, x') = I -J—j-Ax - x'f 1 + 0((x - x'Y) (10) 



f2 1 x+x 
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It is worth point out that the anisotropic polarization factor Eq.(8) still depends on the 
wave field proper-time s. However in this "short-time arrivals" regime, one can expect that 
the proper-time integration can be considered as an average and thus can be factorized in our 



^This WKB's short-time arrival may be thought as the proper-time vanishing hmit 5 — > 0+ on the path 
integral weight: 



exp —77; / d'y 



1 /dr" 



45 7o 'C2(r(7))Vrf7 



4 



general Eq.(4) as a first significant leading approximation for the proper-time integratio: 
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ds Eq.(6) X Eq.(7) 
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In this effective case Eq.(ll), one can easily disregard the "non-physical" advanced piece of 
the soluton as we show in the simple example of a constant velocity medium in R^, procedure 
to be done by hand, since there is an causal ie prescription already bult in. See footnote 3. 
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(12) 

Finally, let us remark that Eq.(4)-Eq.(5) present interesting prospects for implementing 
numerical path-summation methods as exposed fully in refs. [3], see also Appendic C for a 
toy- model of ours to understand the principles of path- integral tomography mechanism ([4]). 



•^This kind of wave propagation factorized form for the wave field at the WKB short-time arrival, where the 
anisotropic polarization vectors dynamics is factorized out from the "scalar propagation ' ' is naturally expected, 
i.e.: 

/>oo 

^,{x)^j{x')^ / ds{5,j +<S>,j){x,x'), 

JQ 

with {^k{x)}k=i,2,3 denoting the wave polarization vectors. 
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APPENDIX A 



Elastic Wave Propagation — Some Comments 



The usual wave equation for modeling an isotropic, non-attenuating, isotropic, linearly whole 
space is given by the complete anisotropic wave equation 

^U%x,t) = A,.,,(x)^-^[/'=(x,i) + DM{^)—U\x,t) + Mx,t) (A-1) 

where the spatially variable anisotropic elastic tensors Aj^j.^ and Dj^^ are explicitly given in 
terms of the medium elastic parameters, namely: the density p{x) and the Lame Parameters 

and X{x) 

A«m(i) = {6"6"'A(x) + {f^){x)(S"S"') + {!^){x)(S"S"-) 

^ ( (*«) + ( -M) (^") + i (A-2) 



p / \p dx^ J \p dx^ 

Let us make a detour on the usual approach and point out that the vectorial variable change 

U\x,t)^[M]k,{x)(t>%x,t) (A-3) 

where the matrix [M]kqi l<A;<3,l<g<3is inversible everywhere and uniquely determined 
from the relationship PDF's first order matrixial system written below 

K„H{x)^^Mks{x) + Kjkr{x)-^Mks{x) + Dikr{x)Mks{x) \ ^{^^t) = 0- (A-4) 

This changes our full elastic wave equation Eq.(A-l) into a more suitable equation of the 
form without the troublesome first-order matrix valued derivative term 

^4>Kx,t) = n,^^,{x)^^^n^,t) + V{x)nx,t) (A-5) 



where the new effective anisotropic tensor is now exphcitly given by 

— 1i / \/ A /\\r / \ (Psrturbative Hypothesis) _(q) , 

and the scalar potential 



X] 



(A-6) 



V{x) = [M-%{x) I (a.,u {g^iM,,^ ) + (D^H^,^'h}j I {x). (A-7) 

One now can proceed as in the bulk of our Brief Report in order to produce a space-time 
path-integral representation for Eq.(A-5). 

As a last comment with a computer oriented focus, some times it appears more suitable to 
write the elastic wave equation into a form of a hyperbolic system involving the propagation of 
a compressional wave 9{x,t) and the shear wave = |^ — namely: we have the set 



p0u = [iAX)e + 2djXdje + XAe + {djdiii)U'^ 

+ {diiJi){die + AC/0 + {djlAi^Uj + djO) + 2AtA^] 
+ 5{x -x)5{t-t) 
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(A-8) 

(A-9) 
(A-10) 



4n\x — x'\ 

puit = ^Jii^u' + xdiO + ^Ji^ie + diXe + ^j^Ji^iU^ + p 

where the source p into Eq.(A-lO) must be compatible with the compressional explosive be- 
havior delta function source in Eq.(A-8). Namely 
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5{t-t'). 



(A-11) 



Numerical Analysis of Eq.(A-8) / Eq.(A-lO) may be implementable by putting the medium 
in a Box and considering truncated Fourier expansion for all spatially variable objects involved 
and leading thus to a set of time-dependent ODE's for the dynamical objects Fourier coeficients. 

For a compact volume Q and the spectral eigenvalue and eigenfunctions of the Laplacian 
Operator in (with Dirichlet conditions) and n e Z"*" fixed 
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(A-12) 
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APPENDIX B 



Feynman Short-time amplitude 



In this appendix we present our proposal proper-time Feynman asymptotic expansion for 
the anisotropic wave operator Eq.(3). We have thus the chains of asymptotic equahties (still 
mathematically formal, unless the operator L is essentially self-adjoint and acting on a spectral 
- finite spectral range parameter - sub-space of £, which means that £ is a bounded operator 
there!): 



lim ((x, t, ^)| exp(zs>C)|(a:', t', A;)) = I{{x^t,i); {x' ,t' ,k)) 




d^pdw{exp{is[-{w'^)6pq + Cpjqi{x)pjPe + w{x)w]})ik 



X exp(ip(a; — x')) x exp{iw{t — t')) 




X{5nk + ise-fnukv{x)PuPv)- 



We have thus (see ref.[l], Chapter 24, Eqs. (24.11) and (24.16)) 



I{{x,t,ty, {x',t',k)) {(C2(x))-^/2(exp ( - ^ 
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Note that we have used the wearly smallness anisotropic perturbation condition to arrive 
at the above written final asymptotic short-times result Eq.(B-l) 



<<£ 



Snk - ^^^,2^'^ ^ )i^uv) ^ ^nk (B-2) 
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APPENDIX C 



A toy model for tomographic seismic inversion in scalar 
path integral wave propagation-highlights 



Let us consider for simplicity of our exposition scalar waves of a definite frequency w (see 
Chapter 20- [1] for the notation details) in a medium with a spatially slowly varying index of 
refraction r]{x) in R^. Namely 

where we have introduced the L. Schulman somewhat artificial scale B > 0. 

Let us further consider the "perturbative ' ' decomposition of the medium material proper- 
ties (with e << 1) 

reference perturbation 

' ^^Z^ (C-2) 



(n(x))^ (no(a;))^ (ni(a:))^ 
where one supposes that all wave dynamics relative to the reference unperturbed medium no(x) 
is known exactly. 

The scattered field by the medium ni{x) is given in the lowest e order by the following path 
integral representation 

U^^{xs,Xr, [n{x)]) 



Jr{0)=xs 

0))) 
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The problem of determining the perturbing refraction index ni{x) - after 
"SampUng" the scattered field U^'-'{xs,Xr, [n]) at {2k + 1) points Xg = 



(X 



(s) 
1 1 



Xr = {xi \ . . . X2,!_^_{), namely: 
One knows the set 



l<g<2A;+l 

at the source Xg*'' and receiver Xp ^ grid point of our Geophysical surveilence, lead us to the 
seismic tomographic scalar inversion problem in the path integral framework in order to estimate 

We, thus, consider firstly the Fourier Transformed of the refraction perturbing index (its 
square inverse for sure!) 

(;^) ^(S^ 

After considering the sampling of V^(^) at N — 2k{k + 1) points; {V{k-N), ■ ■ ■ ■,V{kN)} 



still to be determined from the linear equation (C-3) namely {2k + 1)^ = 
(2 • 2k{k + 1)) + 1), one obtains the short-time "staching' ' equations ([4]). 
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\ TjSL( Jr) is) j 

with the exactly evaluated path integral coeficients sources-receivers propagations as given 
below by Eq.(C-6). (One can evaluate exactly them for instance by means of the short-proper 
time limit 5 — > oo in Schulman's approach, however by keeping limg_»o+ [£( ^.!-f^ )] < -|-oo) 
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Note that the well-posed solution of the linear matrix equation (C-5) is given by the standard 
methods of Linear Algebra 



Finally the problem of handhng Caustics and WkB next corrections with a non-trivial 
determinant associated to the density of paths, can be in principle straightforwarde implemented 
in our proposed tomographic path integral scheme by just considering one- loop corresctions to 
the path integral source- receiver propagator Eq.(C-6), namely, the next one- loop correction 
leading to the Van-Vleck-Morrete determinant is coming from the following quadratic path 
integral 



where the geometrical ray path classical trajectory is given by the well know ray equation 



[V] = {[Mfmr^iMf . [u^^]) 



(C-7) 



After evaluating the FFT coeficients {V{k-j)}N<j<N, one can "estimate' ' Eq.(C-4). 
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